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I I Abstract 

O |. We consider an ensemble of interacting charged particles on the line consisting of 

two species of particles with charge ratio 2 : 1 in the presence of the harmonic oscillator 
potential. The system is assumed to be at temperature corresponding to /? = 1 and the 

C^ ■ sum of the charges is fixed. We investigate the distribution of the number as well as the 

spatial density of each species of particle in the limit as the total charge increases to 
oo. These results will follow from the fact that the system of particles forms a PfafHan 
point process. We produce the skew-orthogonal polynomials necessary to simplify the 
related matrix kernels. 

> 

C<i ' 1 Introduction 

(N 

[^ , Let L, M and N be non-negative integers so that L + 2M = N, and consider l-dimensional 

(^ ' electrostatic system consisting of L particles with unit charge and M particles with charge 2. 

^D . We will identify the state of the system by pairs of finite subsets of M, ^i — {ai,a2, ■ . ■ , ai} 

and ^2 = {/3i, /32, . . . /3a/}, where ai,a2, ■ ■ ■ ,aL represent the locations of the charge 1 
particles and /3i, /32, . • . , /?m represent the locations of the charge 2 particles. 
The potential energy of state ^ — (^1,^2) is given by 



^ 



L M 

^log|a, -afc|-(-4^ log|/3™-/3„|-f 2^^ |a£-/3™|. 

j<k m<n 1=1 m=l 

We assume that the system is in the presence of an external field, so that the interaction 
energy between the charges and the field is given by 

L M 



^F(a,)-2^T/(/3,„) 



for some potential F : R — > [0,oo). Eventually we will specify to the situation where V is 
the harmonic oscillator potential, but for now we maintain generality. The total potential 
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energy of the system is therefore 



M 



E = Y,^og \aj - afel + 4 ^ log |/?„ - /?„! + 2^ ^ log \a, ~ /3„,| (1.1) 

2<k m<n £—1 m—1 



L M 



^yM~2^y(/3„). 



1=1 m=l 

Given a pair of vectors (q:,/3) G M.^ x M^^ we will define E{a.,f3) to be the right hand side 
of (JLII) . and call {a.,/3) a state vector corresponding to the state £,. Generically, there are 
LlMl state vectors corresponding to a given state. 

Assuming the system is placed in a heat bath corresponding to inverse temperature 
parameter /3 = 1, then the Boltzmann factor for the state vector (q:,/3) is given by 

L M L M 

g-£(«,/3)^-Q^(^^) ^w{Prr.Y\{\a,-ak\ J] l/5m-/3n|'nn \'^i - Prn\\ (1.2) 
1=1 m—1 j<k m<n £=1 m — 1 

where ^(7) = e~^^'^' is the weight of the system. The partition function of the system is 
given by 

Zl.m = j^J f e-^(«-^)dA.^(a)d/^(/3), (1.3) 

where /z and ^^ are Lebesgue measure on R and M^ respectively. The multiplicative pref- 
actor 1/{L\M\) compensates for the multitude of state vectors associated to each state. 

Here we will be interested in a form of the grand canonical ensemble conditioned so that 
the sum of the charges equals N. That is, we consider the union of all two component 
ensembles with L particles of charge 1 and M particles of charge 2 over all pairs of non- 
negative integers L and M with L + 2M ~ N. The partition function of this ensemble is 
given by 

ZiX)= Y: X^Z,,m= E^ [ [ e-^(-/^)d^^(a)d/^(/3). 

(L,M) (.L.AI) * " 

Here X > is the fugacity of the system, a parameter which controls the probability that 
the system has a particular population vector {L,M). The sum over {L,M) indicates that 
we are summing over all pairs of non-negative integers such that L + 2M = N. 

Note now that {L,M) is itself a random vector, though we will continue to use this 
notation for the value of the population vector as well. For example, for each admissible 
pair {L,M), the joint density of particles given population vector {L,M) is given by the 
normalized Boltzmann factor, 

^^e-^t"^'^). (1.4) 

Z{X) ^ ' 

When X = 1 the probability of seeing a particular pair [L, M), or Prob(i, M), is the ratio 
Zl,m/Z, where Z = Z{1). 

Experts of random matrix theory will have already noticed that when X — the above 
reduces to a general orthogonal (or (3 — 1) ensemble. Likewise, as X — )■ 00, the above 
formally goes over to the corresponding symplectic (or /3 = 4) ensemble. This provides then 
an unusual sort of interpolation between two classical and well studied point processes. 
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2 Statement of results 

In this paper we will primarily be concerned with global statistics of the particles when the 
fugacity equals 1 and the potential V is given by 



Vij) = 7V2, 



that is, 



w{j) 



-7V2 



Many of the results presented here are valid for other potentials and other values of X, 
however unless otherwise indicated we will restrict ourselves to these choices of V and X. 
We will also restrict ourselves to the situation where N = 2Jis an even integer. 

Similar results for the two-charge ensemble constrained to the circle with uniform weight 
were obtained by P.J. Forrester (see 5.9 of [7] and the references therein). 

The goal of this paper is to present global results about the distribution of L and M 
as well as the global spatial distribution of each of the species of particles. Along the way 
we will derive a PfafHan point process for the particles (similar to that of another two- 
component ensemble, Ginibre's real ensemble) as well as the skew-orthogonal polynomials 
which allow us to present a simplified matrix kernel for the process. The local analysis of 
this kernel {i.e. its scaling limits in the bulk and at the edge) as well an investigation of the 
right-most particle of each species will appear in a forthcoming publication. 

2.1 Distribution of the population vectors 

Sharp results on the law of the state vector (L, M) are consequences of the following char- 
acterization. 

Theorem 2.1. For each non-negative integer j, let Lj = Lj he the jth Laguerre 

polynomial with parameter a — —1/2. Then, PToh{L, M) is the coefficient of X^ of the 
polynomial iyjv/2(~^^)/-^Af/2(— !)• That is. 



Z{X) _ LjV/2(-^') 

Z Li^/2{-l) 



and so 



1. Prob(i,Af) = 



E 



L\M\ '^ l\m\ 



if L is even, and is equal to otherwise. 



2. E{V'' 



^ dx' Ljv/2(-l) 



for m non-negative integer. 



x=i 



Properties of the Laguerre polynomials now allow for nice expressions for the mean, 
variance, etc. of L for all finite values of N. For example, we have that 



•^[^l - dX 



Z{X) 



z 



x=i 



N/2-1 I- 



jv/2-1 



2E 



Lj{-1) 



2E 



E 



i=0 



r 



N 



i\V[i + 



N 



-i+l ]T\i 



Asymptotic descriptions of the law of L are just as readily obtained from Theorem 12.1 
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Theorem 2.2. As N ^ oo it holds: 

1. E{L) ^ V2N - 1 + -^ + 0{N-^) and Var(i) ^ V2N - | + 0{N-^/^), 



3VN 

L- (2iV)i/2 

2. — :-rj2 — converges in distribution to a standard Normal random variable, 

3. Prob (|^= - 1| > e) < CNe^^''^'^''^^^ with a numerical constant C for any e > 0. 

2.2 Spatial density of particles 

We introduce the (mean) counting measures pi and p2 for the charge 1 and charge 2 particles 
defined by 

^[\Ar\ii\]^ ( dpi and E[|An6|] = / dp2 

J A J A 

for Borel subsets A C M (where, for instance, |A n fi| is the number of charge 1 particles 
in A). As we shall see in the sequel, these measures are absolutely continuous with respect 
to Lebesgue measure, and we will write R\ q {x) and Rq i [x) for their respective densities. 



(The cryptic notation will be resolved in Section 13. 2[ when we define the t, rn-correlation 

function of the ensemble to be R^J^)- 

From Theorem l2.2l we see that, as A^ — > oo, 

r[%^ (x) dx = E[L] - V2N, and / i?[,^^ (x) dx = E[A/] - — . 

-OO J —OO ^ 

One then would ask, when suitably scaled and normalized as in 



(x) = ^i?i7(VAx) and sr^(x) = — i?^;V(VAx), 

whether s]^ (x)dx and S2 (x)(ix converge to proper probability measures. This is answered 
in the affirmative in Theorem 12.31 below. 

The previous result shows that, with probability one, for all N large the number of charge 
1 particles is v2]V(l+o(l)). This suggests that, in the thermodynamic limit, the statistics of 
the charge 2 particles should behave as though there are no charge 1 particles present, or like 
a copy of the Gaussian Symplectic Ensemble (again, arrived at from the present ensemble 
upon setting L = 0). Indeed we find the scaled density of charge 2 particles approaches the 
semi-circle law. 

On the other hand, though the charge 1 particles exhibit the same level repulsion amongst 
themselves as the eigenvalues in the Gaussian Orthogonal Ensemble (occurring here when 
M = 0), the preponderance of charge 2 particles leads to a different limit distribution. 

Ginibre's real ensemble, the ensemble of eigenvalues of real asymmetric matrices with 
i.i.d. Gaussian entries, has superficial resemblance to the ensemble we are considering here. 
First, it is suggestive to think of the present ensemble as arising from real Ginibre by forcing 
the non-real eigenvalues, which occur in complex conjugate pairs, to be identified with one 
"charge two" particle on the line. A little more concretely, the (random) number of real 
eigenvalues in real Ginibre has both expectation and variance of 0(vA), as does the number 
of charge 1 particles here. (See [6] for the mean, and [8] for the variance). It is perhaps not 
surprising, therefore, that the limiting scaled density of charge 1 particles is the same (up 
to a constant) as that of the real eigenvalues in Ginibre's real ensemble [4]. 
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Theorem 2.3. As N —>■ cx), s\ converges weakly in the sense of measures to the uniform 

law on [— v2, v2]; cind Sj converges in the same manner to the semi-circular law with the 
same support. In particular, it is proved that 



Jtx (N) 



{x)dx 



V2t 



sm{V2t) 



Jtx(N) 



{x)dx 



^Ji(V2t), 



where the convergence is pointwise. 

We give an elementary proof of the above, making use of the exphcit skew-orthogonal 
polynomial system derived below. Given that the number of charge 1 particles is o(N) , one 
could undoubtedly make a large deviation proof along the lines of [2] or [5] of a stronger 
version of the second statement: that the random counting measure of charge 2 particles 
converges almost surely to the semi-circle law. However, it is not clear how to use such 
energy optimization ideas to access the charge 1 profile. 





„W 



„(^ 



Figure 1: s]^ (left) and 53 (right) for, from lightest to darkest, N — 10,30 and 90. 



3 A PfafRan point process for the particles 

All of the results in this paper follow, in one way or another, from the fact that our interacting 
particles form a Pfafhan point process very much like that of Ginibre's real ensemble and 
related to the Gaussian Orthogonal and Symplectic Ensembles. 

The results in this section are valid for quite general weight functions w and fugacities. 
Thus, for the time being, we will return to the general situation. 



3.1 The joint density of particles 

The joint density of particles for a particular choice of {L, M) is given by 
X^ 



Z{X) 



nLMia,(3), where f]L,M(a,/3) = e"^'"-'^). 



More specifically. 



M 



L M 



nL,Mia,(3)^Y[wiae) ]]_ wiP^f Ula^ - a^] J] Wm-pnl^H J] 1"^ " /^^ 



m\ ) 



e=i 



J<k 



£=1 m=l 



where, for now, the only assumptions we will make on w are that it is positive and Lebesgue 
measurable with < Z{X) < 00. 
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3.2 Correlation Functions 

Given < £ < L and < m < M, wc define the i, m-correlation function i?^ J : K^ x M™ 

[0,(X))by 



<^(x;y)= E 



{L,M) 
L>£,M>-n 



{L - £)\{M - m)\ 



f]L,M(xVa,yV/3) dAi'^-^Ca) ^''"'"(/S), 



wfiere, for instance, xVa is tlie vector in M-^ formed by concatenating x e M^ and a E R^ ^. 
We will often write Re^m for R\ J in situations where N is seen as being fixed. 

The correlation functions encode statistical information about the configurations of the 
charged particles. To be more precise, given a E M^ and /3 e K^^ with L + 2M = N, we set 

e = e(a, f3) - (Ci, 6) - (?i(a), ?2(/3)) = ({«!, • • • , aj, {/3i, . . . /3™}) . 

Given an i-tuple of mutually disjoint subsets of R, A = {Ai, A2, . . . , Al), and an M-tuple 
of mutually disjoint subsets of R, B = (i?i, B2, . . . , Bm), the probability that the system 
is in a state where there is exactly one charge 1 particle in each of the Ai and exactly one 
charge 2 particle in each of the Bm is given by 

Prob{|Ai n Cil = 1, . . . , |Al n 61 = 1, |Si n 61 = 1, • • • , I^M n 61 = 1} 



= E 



1=1 



\{\Atr\iA\ n |5™nC: 



M 



This probability can also be represented by 
TaJi ^ ^ 



L\M\ 



Oi,M(a,/3)rfM'^(a)V(/3) 



Since the integrand is symmetric in the coordinates of a and /3, we find 

.^1=1 ^^m=l J J -^B^A 



OVn/3) 



The correlation functions can be used to generalize this formula. If A = (Ai , yl2 , . . . , A^) 
is a tuple of disjoint subsets of R and B — {Bi, B2, . . . , Bm) another such tuple, then 



E 



j=i ^ *- fc=i 



i?,,,„(x;y)V(a)ci/"(/3). 



B J A 



3.3 PfafRan point processes 

Consider, for the moment, a simplified system of indistinguishable random points Q 

{71, 72, ... , 7Ar} C R with correlation functions Rniz) satisfying 



E 






Al 



i?„(z)d/i"(z) 
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for any n-tuple {Ai, A2, . . . , A„) of mutually disjoint sets. 

If there exists a matrix valued function Km -.IS.^ ^^^^^"^ such that 



then we say that our ensemble of random points forms a Pfaffian point process with matrix 
kernel K^. Much of the information about probabilities of locations of particles {e.g. gap 
probabilities) can be derived from properties of the matrix kernel. Moreover, in many 
instances, we are interested in statistical properties of the particles as their number (or 
some related parameter) tends toward 00. In these instances, it is sometimes possible to 
analyze Kpf{x,y) in this limit (under, perhaps, some scaling of x and y dependent on N) 
so that the relevant limiting probabilities are attainable from this limiting kernel. 

For the ensemble of charge 1 and charge 2 particles with total charge N, we will demon- 
strate that the correlation functions have a Pfaffian formulation of the form, 



i?£,m(x;y)=2^Pf 



K 



jY \Xj,Xj') ^ j\f \XjiZk') 

K'Jl^{zk,Xk') K'jf{xk,Xk') 



k,k' 



1,2, 
1,2, 



.m 



where Kj^ , Kj^ , Kj^ and Kj^ are 2x2 matrix kernels. 

3.4 A PfafRan form for the total partition function 

In order to establish the existence of the matrix kernels we first need a Pfaffian formulation 
of the total partition function. 

Given a measure j^ on M we define the operators e^ and 63 on L'^{v) by 



ei/(a;) = 2 / /(y) sgn(2/ - a;)dzy(2/) 



and 



£2/(2/) - /'(y). 



(Obviously £2 does not depend on v, but it is convenient to maintain symmetric notation). 
Using these inner products we define 



(/l5) 



62 



[.f{x)ebg{x) - g{x)ebf{x)] dv{x), 



1,2. 



We specialize these operators and inner products for Lebesgue measure /i by setting e^ ~ e^ 
and {f\g)'^2- We also write f{x) — w{x)f{x). It is easily seen that 



(/l5)i == / f{x)eig{x) - g{x)eif{x) dfi{x) = {f\g)l 



•.wn 



Similarly, 



(./l5)4 = 



f{x)^gix) ^ gix)^f{x) 



dfi{x) 



wixf [/(x)g'(x) - g(x)/'(x)] dx ^ {f\g) 



W fl 

4 



We call a family of polynomials, p = {pq{x),pi{x), . . . ,pM~i{x)j , a complete family of 
polynomials if degp„ — n. K complete family of monic polynomials is defined accordingly. 
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Theorem 3.1. Suppose N is even and p is any complete family of monic polynomials. 
Then, 

Z(X) = Pf (X^AP + BP), 

where 



AP 



lAf-l 
^in\Pn/l\,n.7i=0 



i\Pn 



and BP 



i\P' 



iN-l 



■m\Pn/4:\jn,i=0 ■ 



Corollary 3.2. With the same assumptions as Theorem \3.1\ Z = Pf(AP + BP). 

3.5 A PfafRan formulation of the correlation functions 

In order to describe the entries in the kernels Kj^ , Kj^ , Kj^ and Kj^ , we suppose p is 
any complete family of polynomials and define 

CP = AP + BP, 



where Ap and Bp are as in Corollary 13.21 (Here we are setting X — 1, though similar 
maneuvers are valid for general X > 0). Since we are assuming that Z = Pf Cp is non-zero, 
CP is invertible and we set 

The Cj.fc clearly depend on our choice of polynomials. We then define 



JV-l 

>CN{x,y) = ^ Pj{x)Q^kPk{y)- 

j,k=0 



(3.1) 



The operators ei and £2 operate on >CN{x,y) in the usual manner. For instance, 



and 



JV-l 

e2XN{x,y) ^ ^ e2P3{x)CjM£2Pk{y) 

j,k=0 

XNei{x,y) = ^ Pj{x)Q,keiPk{y)- 



(That is, e written on the left acts on the >CN{x,y) viewed as a function of x, etc.). 

Theorem 3.3. Suppose N is even, p is any complete family of polynomials and >CN(x,y) 
is given as in I13.1\) . Then, 



Re,mi^;y)^2'Pi 



K'^^{zk,Xk') K'^'^{xk,Xk') 



jj' = l,2,...,£ 
k,k' = 1, 2, ... to' 



where 



Kpi i^^y) 



>CN{x,y) XNei{x,y) 

ei>CN{x,y) eiKNei{x,y) + jSgn{y-x) 



K'^/{x,y) = 



XN{x,y) XNe2{x,y) 
e2>CN{x,y) e2>CNe2{x,y) 



K]f{x,y) 



XN{x,y) KNei{x,y) 
e2XN{x,y) e2XNei{x,y) 



and K'^^{x,y) 



>CN{x,y) >CNe2(x,y) 
eiXN{x,y) ei>CNe2{x,y) 
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Remark. The factor 2^ can be moved inside the PfafBan so that the entries in the various 
kernels where an ei appears are muhiphed by 2. This maneuver is superficial, but has 
the effect of making these particular entries appear more like the entries in other /3 — 1 
ensembles (e.g. GOE). For instance, ^ sgn(j/ — x) appears more natural to experts used to 
these other ensembles. 

We can simplify the presentation of the matrix kernels with a bit of notation. First, let 
us write 



KN(x,y) 



}<N{x,y) >CN{x,y) 
XN{x,y) XN{x,y) 



and 



Eb 



l 0' 
efc 



1.2. 



Then, 



'1,1 



Kj^ ix,y) ^ EiKN{x,y)Ei 











i sgn{y ~ x) 



KN^{x,y) = E2KN{x,y)E2, K]f'{x,y) == EiKN{x,y)E2, K'^^^{x,y) = E2KN{x,y)Ei. 



by 



We notice in particular that the functions R[ q and Rq ^ given in Section [2^ are given 



N-l 



Af-1 



j.k=0 



and Ro!i{x) = ^ Pj{x)Cj.k£2Pk{x). (3.2) 

3,k=Q 



3.6 Skew-orthogonal polynomials 

The entries in the kernel themselves can be simplified (or at least presented in a simplified 
form) by a judicious choice of p. If we define 



(fig) - (/l.9)i + (/l5)4, 



then 



CP = [(p™|p„ 



Jr?i.n— ' 



Since ^c^ (and by extension all other entries of the various kernels) depend on the inverse 
transpose of C^, it is desirable to find a complete family of polynomials for which C^ can 
be easily inverted. 

We say p = {po,pi, . . .) is a family of skew- orthogonal polynomials for the skew-inner 
product ('I-) with weight w if there exists real numbers (called normalizations) ri,r2,... 
such that 

{p2j\P2k) = {p2j+l\P2k+l) = and (p2]\P2k+l) = -{p2k+l\P2]) = Sj,krj. 

Using these polynomials, the entries in the matrix kernels presented in Section [3.51 have 
a particularly simple form. For instance. 



j-i 

XN{x,y) = Y^ 



P2j {x)p2j+i [y) - P2J+1 {x)p2j (y) 



and the entries of the kernels are computed by applying the appropriate e operators to this 
expression. 
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3.7 Specification to the Harmonic Oscillator Potential 

We now return to the case where the weight function is w{x) = e^^ '^. 
Theorem 3.4. Let 

(.|.)m=X2(.|.)l + (.|.)4. 

A complete family of skew- orthogonal polynom,ials for the weight w with respect to (•|-)'-'''^-' 
is given by 

^if(-) = E(-l)'^^T7^^^(-')' (3.3) 



^^(0) 



and 



p!^^,{x) ^ 2xPif\x) - 2^P^f\x) 

= 4X^xVi-lf ^^'^7^'h lix^)+2xi-ir ^"'^^7^'h iix'). (3.4) 

fe=o LliO) L„^{0) 

where Lk{x) — L\. (x) is the generalized kth Laguerre polynomial. The normalization of 

this family of polynomials is given by 

(H^^'lP^lr^ = '^^^L^i-X^)L^M-X^ (3.5) 

r(m+ 2) 

We can recover a family of monic skew-orthogonal polynomials by dividing by the leading 
coefficient. Specifically, 

Corollary 3.5. A complete family of monic skew- orthogonal polynomials for the weight w 
with respect to (•l-)'-'''^-' is given by 



an- 



d 



(X) f s iX)f y "^ (X)/ \ 

The normalization for this family of monic skew- orthogonal polynomials is given by 
^ix)_^x)^x) gx) (.7 + l)!r(7 + |)L,+,(-X^) 

Setting X = I, we recover a family of skew-orthogonal polynomials for the harmonic 

oscillator two charge ensemble with fugacity equal to one, and we will write p„ for p„ and 

f (1) 
rj tor r- . 
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4 Proofs 

4.1 Proof of Theorem \2A] 

We set J = N/2. To prove 1, we use Theorem 13.11 and the skew-orthogonal polynomials 
from Corollarv l3.5l to write 



Z{X) = Pf 



r, 



{X) 



r 



J-i 



J-i 

J=0 



„(^) 



Hence, 



Z{X) -'^rf^ Lj{~X^)Lo{-l) Lj{-X^) 



n 



_ ^0 r, L,(-l)Lo(-X2) LA-l) ' 

where again Lj{x) — Lj ' {x). Note Lo{X) = 1. 

The remaining claims follow from the above by definition and the properties of Laguerre 
polynomials. 



4.2 Proof of Theorem [272] 

Point 3 of Theorem 12.11 specified to the first two moments produces 



E\L] = 



dX 



Z{X) 

z 



Var(i) = 



x=i 

1/2 



dX\ dX Z ) \dX z ) 

1/2/ N ri/2 



J x=i 



Now, since L'j{x) — ~L j_-^{x) and Lj{x) — Lj{x) — Lj_^{x), we have that 

E(L) = 2 



.4/-l(-l)_o4''(-l) 



L,,(-l) Lji-l) 

Further, using the differential equation xL'j{x) + (1/2 — x)L'j{x) + JLj{x) — 0, we also have 
that 



d 
dx 



(x^Lj(-x^)\ = -4xL'(-x^)+4x^L"j{-x^) 
\ dx J 

= {-2x + 4:X^)L'j{-x^) + 4:xJLj{-x^). 



This yields 



Var(L) = 4J - E{L) - E{Ly 



and so asymptotics of the variance follow from those for the mean. 
Next introduce a version of Perron's formula (see [5]), 

L^(-l) = -^m"/2-i/4g2v^ /^ ^ Ci(a)m-i/2 + ^2(0)771-^ + 0(m-=^/2)) 
2v7re V / 
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where m = n + 1 and Cj{a) are known explicitly. In particular, Ci(l/2) ~ —1/6, C2(l/2) = 
-7/144, Ci(-l/2) = -2/3, and C2(-l/2) = 77/144. Substituting into the above we then 
obtain 

E(L) = 2V7TT - 1 - g^l^ + O(J-i) = 2x/j - 1 + ^ + O(J-i), 

and Var(L) = 2\/j — | + 0( J^^^^) which completes the proof of point 1 (recall J — N/2). 
Moving to the limit law for L, we introduce a little new notation. Set 

^"^'^ ^ r(f-| + i)^ + i)r(| + i) = ^"'^"(')" 

with Cat = r(Y + 5)[-^Af/2(— 1)]^"^- For A; even, pAr(fc) is the probability of k particles 
of charge 1, otherwise this probability is zero, compare point 1 of Theorem 12.11 In the 
continuum limit this distinction is unimportant; we will show that, as A'^ — > oo 



— c /2 

(2iV)i/4pjv((2iV)i/2 + (2iV)i/4c) = ^^(1 + OiN-^/^)) (4.1) 



/2^ 
uniformly for c on compact sets. 

First note that by Stirling's approximation (in the form r(z) = w ^(z/e)^(l + 0{j))) 
and again Perron's formula (now in the simpler form Lz{—1) — 2 /Irez ^"^ ^^ '^ ^(~7?)))' 

Cn = 2ttVn^ (7V/2) W2)g-JV/2-y2]V(^ _^ 0(7V-i/2)). (4.2) 

Next, with both k and N — k large we have 

QNik) = (2^)3/2 VM(7V/2)(^/2)g-Ar/2-fc/2 (43) 

X g[(iV/2-fc/2)log(l-fe/iV) + (fc/2)log(feV2W)l(l + 0(^-1 V (TV _ fc)-l V fciV-1)), 

again by Stirling's approximation. Restricting to fc = 0{^/N), (|4.2p and (|4.3p yield 

^^(^^ ^ y^e-*"^'Hl + 0(7V-i/2)), (4.4) 

where 

Mfc) = ^/2iV-- + (---)log(l--) + -log(-). 

Now, quite simply 

) log(l = + + 0(N-^^^), 

2 2-^ ^^ n' 2 an ^ '' 

if fc = 0(\/]V), and, if k is also such that 1 — 557 ~ 0(A^^^/^), we further have 

1 log(— ) = --(1 - — ) - -(1 - — )2 + 0(iV-i/4). 

2 ^4iV^ 2^ 2iV^ 4^ 2N' ^ ' 

More precisely, from the last two displays we readily find that 

I e2 



•n{V2N + i) = ^+ " + 0{N-^/-^), uniformly for i = 0{N^^'^) 
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Substituting back into dO]), since {V2N + i^^/^ ^ (2iV)-i/4(l + 0(Ar-i/4)) ^g^^-^^ f^^. 
£ = 0(Ar-V4)^ completes the verification of (|4T1) . 

Last, for the tail estimate, revisiting (|4.2p and (|4.3p shows the conclusion of (|4.4p may 
be modified to read 

for all 1 < fc < A^ with a numerical constant C. (Here we understand (1 — ■^) log(l — ■^) to 
be zero at k — N.) Differentiating yields 

A^^(fc) = ilog(^^^^^ 

and so (j)N{k) is decreasing for k < c^^^/2n and increasing for k > c\/2N for any c > 1. 
Now, since (1 - e) log(l - e) > -e and log(l + e) > e - e^/2 for < e < 1, 

0Ar((l + e)\/2]V) > -e\/2iV + 2(l + e)\/2iVlog(l + e) >e\/2]V, 

also for < e < 1. Hence, for c > 1, Prob(L > cV^) < NpNicV^N) < 
CNe-^^"-^^^^^"^^. The proof for the left tail is much the same. 

4.3 Proof of Theorem [2731 

In both cases we use the expression of the one point function in terms of Hermite polyno- 
mials, see (4.17) and (4.20) below. 
We start with 

(1). N _ /T^ V^ eiP2n+lix)p2n{VNx) - P2n+l{VNx)eiP2n{VNx) 

^N (X) — V-i 2_^ , 

n n 

n—0 

and 

2 T^' p'2,,+ Ay^x)P2n{VNx) - P2n+l{VNx)p'2n{^x) 



(2) 



w-T^i: 



N ^ Tn 

n—0 



along with the relations J^^P2n-\-i = eip2n-\-i{x) = 2p2n{x) and ^2ni^) = ^2P2n{x) = 
— ^P2n+i(^)- An integration by parts in both instances then allows: with t^ — t/yN^ 

e'*'' s^^\x)dx (4.5) 

= ^"e'^"' r e^*"^(?2„(x))^dx- ^^''£'r-i r e^'"^p2n{x)e,p2n{x)dx, 

and 

e**^s^^(a;)da; (4.6) 



2 



JV/2-1 ^oo o.. ^/2-l 



E r-^ r e^'"''iP2n+i{x)fdx-^f- ^ r-' r e^'-^p2n{x)p2n+i{x)dx. 

n -'— CO n <y —OQ 



,1 = •'-oo „^o 



The first, and primary, step is to show that the advertised limits stem from the first sums 
on the right of the above expressions. 
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Lemma 4.1. Let s]^ (i) and s]^ (i) denote, respectively, the first term on the right hand 
side of (R3p and ( |^.6'[ ). Then, 



as N ^ oo. 



Proof. Since Hk{x) = Hk{x)e ^ 1"^ are the eigenfunctions of the Fourier transform — in 
particular (ilfc) {x) = — == /^ e^^'^Hk{u)du = i'^Hk{x) — we have that 



P^n{x) = Y.^-lfakH2k{x), (4.7) 

fe=0 

n n 

p^i{x) = i'^{-lfak{H2k+i{x)+^kH2k-i{x))^2ix^(-l)''akH2k{x). 

k=0 fe=0 

The last equality makes use of the three term recurrence Hn+i{x) — 2xHn{x) — 2niJ„_i(a;). 
Plancheral's identity then yields, 



n=0 0<fe,f<n 






and 

8 ^^^"^ r°° 

■^N^t) ^ j;^ J2 "^n' J2 {-^f^'akai {x^-{tN/2f)H2k{x+tN/2)H2t{x-tN/2)dx. 

n=0 0<k,e<n J -<x> 

(4.9) 
We begin with the asymptotic considerations of (|4.8p which is slightly simpler. 
From the expansion Hn{a + 6) = X]fe=o {T)^k(o-){2b)"^'' wc find that 

i/2fe(a; + t/2)H2£(x-i/2)dx = e'^'/^V^ V P^) PV!(-2)"i2'c+2£-2™^ 

m=0 \ / V / 

Given this, s)^ is equivalent, as A^ — ?► oo, to 

n=0 fe=0 m=0 *- ■' V / 

/9^' ^/2-l 2fe . . / +2 \ «-fe+™ 

^27r Y^ -1 Y^ ofe+^/o;MY^ (2l)2l^2k+m / I 



V -^ * -VI— n I — n -m— n \ / \ / 



^ E -n' E «.«.2'=-^^(2fc)! ^ 



^ TO!(2^-2/fc + m)! V 2N 

in which we have introduced a self-evident notation for the diagonal and off-diagonal com- 
ponents as well as the (nontraditional) shorthand (n)m := /„ "^y ■ 

Next, recall the definitions a„ = T^^^n (~1), ^n — \/7r2^"+^(2n + 2)!a„a„+i and 
note the simple appraisals: with c — (47re)~^, 

r-i = ^_e-4V^(l + 0(n-i/2)), (4.II) 

a„a„+™2'"+2™(2n)! = cV^n-"-i/2g4y?r(^ ^ ^(^^-1/2-))^ 
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where the latter will be used for m nonnegative and moderate (compared with rt^'^). We 
will also make repeated use of the fact 



n 

^n"-i/2g4v7l ^ ^mg4v7l(i _^ 0(n-i/2)), (4.12) 



fe=i 



valid for any real m. 

Continuing, we change the order of summation to write 



V^^ m=0 ^ ■' Ti=[m/2] k=\m/2'\ 



Then, for fixed m, 

^/2-l n Afin+1/2 

n—\mJ2\ fc— 1 771/2 1 * v / 

by (|4.1ip and (|4.12p . and a dominated convergence argument yields 



°° f— t^/2y" /■"'" 

lim s^}.>, = V , ,,„,/ ^ ^, = / Jo(V2te)da;, 



(4.13) 



for the diagonal contribution. Next, for the off-diagonal terms (second line of (|4.10p ). we 
again change the order of summation and have s]^^ equal to 

^ g=l m=0 ^ ^ ^ 71=9+^777/2] fc=r777/2] 

With now q and m fixed, 

iV/2-1 n-q 

J2 r-^ Y. (2A: + 2g)2,+7„afeafc+5(2fc)!22'=+9 

n=q+[m/2] fc= [777/2] 

= ^ V ^e-4^ V fc^+"-i/2e4^(l + 0(1)) = ^4 ;(! + 0(1))' 

4V^7l^iV^ ^ 4V2^(2g + 2m+l)' 

again by (|4.1ip and (|4.12p . Hence, for bounded t, 

N/2-1 N-2q-2 (_^2/2)9+m 

^N,o = 2 ^ E rn!(2q + m)!(2g + 2m+l)^^ + °^^^^ 

q—l 771—0 ^ / V -I / 

^ (-^72)^ / 2^^ ^\ 

^ (2£ + i) \{2£y. ii\)y (1 + °^!))' 
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after changing variables and the order ol summation in line two. That is, s]^^ tends to 
^^ — /„ Jf){\/2tx)dx, which, combined with (|4.13p . proves the first statement of the 



72* Jo 

lemma. 

Turning to (|4.9I) . the preceding shows that, asymptotically, the (x^ — (ijv/2)^) within the 
integrand may be replaced by jH2{x) — x^ — 1/2 for which there is the related evaluation: 
assuming k < (., 

H2{x)H2k{x + t/2)H2i{x - t/2)e"'' dx (4.14) 

_ . r- V- /2fc\/2A n!m!2'^t^''+^^-"-"' 

oi^2fc U; W^" ^ (ii_Hi + i)!(iii_ii + i)!(ii±iii_i)!- 

n-m=0,±2 

The resulting diagonal term (when k — t \n ()4.9p ) then reads 

n=l fe=l m=0 ^ '-^ ^ ^ 

TV Z^ » Z^ fc ^ '^ Z^ to! m + 2 ! l '^N 

n=l fe=l m=0 ^ ' / \ 



This object does not converge on its own; cancellations from the off-diagonals are required. 

-(2) 
^p>l ^N,{o,+p) 



f2\ ^(2^ 

With similar notation to the above we decompose s]y ^ as in X]n>i ^n (o + ) '^^ which 



■^N (o +p) ^^ arrived at by choosing ^ = fc + p in (14.91) . Writing out the p — I case in full we 
have that 

.Af/2-l „-l 



(2) ^ _16V^ Y^ „-l\:^„ „ 02/C/0LM V^ (2^ + 2)m+2 /^ t^ 

n—1 k—1 m—0 



.,(o,., - ^^E^«^E«^«-^2-(2fc)!5:if^^r^ -i^) (4.16) 



320F"f-^ "^^ ,2...,.,v (2fe + 2)^+2(2fc-m) / t^ ^ ^^^ 
+ ^r E 'n }_^a,a,+i2 (2fc)!}^ TO!(m + 2)! [-2N 

n=l k=l m=0 \ ' / \ 

16v^ "f-^ ., "- ^-^ (2fe + 2)^,, / t^ X "+^ 

-n^ E '^n ):^a.a.+i2 (2fc)! )^ -^^-^^ (^- — 

n=l fc=l m=0 \ ' ' \ 

Consider now the first sum on the right of (I4.15P for fc = n only: 

Sy/jT Y^ _i 2o2n/o mY^/oM fr, (-i^/2iV)™ 

2^ r„ a„2 (2n)! 2^(2fc)„(2r^-m)- 



n—1 m—0 

- ^E^^ "e" ^-«^2-(2n)!(2n)™+^(l + o(l)) 

m=0 ^ ' n=[m/2] 

- 2f:-^^^g^(i+o(i)), 

"^-^ (m!)^(m + 1) 

m—0 ^ ' ^ ' 
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by the same type of estimates used in the analysis of s]^ . Next, using the additional fact 
that 1 — 4fcafe_|_iaj7^ = —k~^^'^{l + 0{k~^)) the remainder (or fc < n — 1 part) of the first 
sum in (|4.15p plus the first sum in (|4.16p is asymptotic to 



8^^^'i-ty2r ^ 



J, E^l^ E r-^ E a^2-(2fc)! X 1 (2.)-^ 

m=0 ^ '^ n=|"m/2] fe=[m/2] ^ 

The last two displays combine to produce the advertised limit ^Ji{\^t). 

^(2) ^(2) 

The above ideas propagate. In particular, the remaining terms of s)^ d~^^N (o +i) balance 
to produce a o(l) contribution, and this appraisal extends to the full sum over p > 1 of 

(2) 

^N (o +p)- ^^ '^'~' '^°^ reproduce the details. D 

Revisiting second terms in (|4.5p and (14. 6p shows that the proof of Theorem 12.31 can be 
completed by the following (rough) overestimates. 

Lemma 4.2. As N ^ oo, 

Y, r-^ \p2n{x)eiP2n{x)\dx^O{N^/^), Y. ^TW \p2n{x)p2n+l{x)\dx^O{N). 

Proof. Along with the well known evaluation ||i/fe||i2 = 7r^/''2"/^\/nI used several time 
already, it holds that HiJ^H^i — c2"^^\/n\n^^^^{l + 0{n^^) with a (known) numerical 
constant c. Next note that ||/eig||Li < II/IIli||5||li and so 



hneiP2n\\L^ < \\P2n\\'ii = (^ Ofc | |^2fc | |li ) 

fc=0 
n 

< C(^ ^e2v^)2 = 0(6^^) 



k=Q 



Vk 



Here we have used, again, that ak2^^/2k\ ^ fc^^/^'e^^. Recalling that r~^ ^ rT^I'^e^^ 
finishes the first part. 

For the second estimate we can get by with an application of Schwarz's inequality. Simply 
compute 

\\P2n\\h -T.-l\\H2k\\h < ^E 4r^'^ = 0(e^^), 
fc=o k=a V K 

and 

n n—1 

\\p2n+l\\l2 = 5]4(||i/2fc+l||i2+16fc2||iJ2fe_i||2,) + 2^afcafc_l||iJ2fc+l||i2 

fc=0 k=0 

n 

< cY,Vke-'^ = 0{ne^^), 

k=0 

to find that ?'^^||p2n||L2||p2n+i||L2 = 0(1), which suffices. D 
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4.4 Proof of Theorem 13.11 

We will prove something slightly more general which will be useful in the sequel. 
Given measures vi and 1^2 on R, define 






and 






(L,M) 



L,M ■ 



(4.17) 
(4.18) 



Theorem 4.3. Suppose N is even and p is any complete family of monic polynomials. 

Z''''''^{X)=Pi{X^AP + BP), (4.19) 

where 



AP - [{p„.\Pn)?t-l 



and BP = [{PrnlPn)'^] 



N-1 
rn,n—0 ' 



AAA The Confluent Vandermonde Determinant 

A special case of the confluent Vandermonde determinant identity has that 



dot 



1 

a\ ■ 


1 

OLL 

■ «i 


1 

PI 





1 

2/?i 


1 

■ Ph 





1 

2/3m 


.<- . 


• "r^ 


f^N-l 


(N 


-l)/?f-2 • 


• PZ-' 


{N 


" i)/3fr^ 



L M 

= n("fc - "j-) n (/5™ - /3")' n n ("^ - z^™)'- (4.20) 

j<fc m<n 



e=i m=l 



We will denote the matrix on the left hand side of (j4.20p by V(q:,/3) and its determinant 
by A(q:,/3). We will later use the fact that the monomials which appear in the definition 
of V(q;, /3) can be replaced by any family of monic polynomials p — {po,P2, ■ ■ ■ ,Pn-i) with 
degp„ — n without changing the determinant. We will write the resulting matrix VP(a, /3), 
and we note that A(a,/3) = detVP(a,/3). 
It follows from (Oil and (fOl) that 



I dctYP {a, f3)\diy^ (a) dv^' (/3). 



^'^' " L\M\ 



4.4.2 Notation for minors 

Given a non-negative integer L, we denote the set {1,2,...,L} by L. By convention, if 
L = 0, then L is the empty set. Given a function i : L /* iV we denote by t' the unique 
function N — L /• N_ whose range is disjoint from t. We denote by i the function L /• N 
which is the identity on L. 
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We define sgni as follows: Let ei , 62, . . . , e^ be any particular basis for R^. We then 
specify that 

et(i) A • • • A e^(^ A 64/(1) A • • • A et/(jv_i) = sgni • ei A 62 A • • • A ejv- 

That is sgnt = (—1)'^ where k is the number of transpositions necessary to put the set 

iil),---i{L),i\l),---,i\N-L) 

into order. Clearly, sgni — 1. 

Given an increasing function i : L /^ N_ and a vector a G M^, we define the vector 
at e M^ by 

"i = (ai(i)7at(2),---,at(L)) 

If u : L /* iV is another increasing function, and A — [am.n] &n N x N matrix , we define 
At^u to be the L x L minor of A given by 

4.4.3 The Laplace expansion of the determinant 

Using this notation, the Laplace expansion of the determinant is given by 

detA= 2_] sgnt • det At,i • det Af^i'. 

t:L7"A[ 

In particular, the Laplace expansion of the determinant of Vp(q:,/3) is given by 
VP(a,/3)= Y^ sgnt-detVPj(Q;)-detVt/,i/(/3), 

t:L7"JV 

where the notation reflects the fact that minors of the form Vt^j depend only on ct and 
minors of the form Vf,;' only depend on /3. 

4.4.4 The Total Partition Function 

Using the previous deflnitions, we may write 

|detVP(a,/3)|= J2 sgnt|[]sgnK-a,)|detVP^(a)detVP;,(/3), 



and 






1 

M! 



x^ / detVP.i,(/3)d;.*^(/3). 



We define 



Ai = — \Y\ sgn(Q:fe ~aj)\ det Vf.(Q:) dv^{a), 
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and 



M ! J^M 



SO that 

ZZm^ E sgnt-A^St,. (4.21) 

4.4.5 Simplifying Ai 

If L is even, 

]JsgnK-aj)-PfT(a), 

where T(q:) is the L x L antisymmetric matrix given by 

T(a) = [sgn(afe-aj)]^^^^^^. 

A similar formula is available when L is odd, but since N and therefore L is even, we will 
not need that here. 
It follows that 



Ai = ^l PfT(a)-detVfj(a)di^f(a), 



and 



detVPi(Q:) = E sgncr]^pt(£)_i(Q;CT(£)), 
where Sl is the symmetric group of L elements, so that, 

At = — V / <^pt(£)_i(a^(£)) Ugncr •PfT(a)di^f (a). 

Now, Sl acts on R^ by permuting coordinates — denote the action of cr by ct h- > cr • a. It is 
easy to verify that 

Pf T((T • a) = sgncr • Pf T(a). 

Consequently, 

A = i- V / (frpt(,)_iK(,))lpfT(a-a)di.f(a), 
and by reindcxing the integral by o: h- > a^^ • a, we find that 
^t = T?E / |nPt(£)-iM|pfT(a)di.f(a)= / | []pt(,)_i(a,)| Pf T(a) d^^f (a) 

Next, we write L = 2K and expand Pf T(q:) as a sum over the symmetric group: 

1 ^ 

Pf T(q:) = -^Kj^ E '^g^'^ n sgi^("T(2fe) - ar(2fc-l)), 

' tGSi, fe=l 
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so that 



A = i( , y2 sgnr / <^ T\pi{2k-i)-i{a2k-i)pt{2k)-i{a2k)sgn{ar(2k) - aT{2k-i))> di^^{a) 
^ ^- res, ^K"- U=i J 

tESl 
^ / \ T\ PiOT{2k-l)-l{ctr(2k-l))PiOT(2k)-l{ar(2k))sgn{ar(2k) ~ aTi2k-l))\ dl^tioi) 

1 ^ /• 

" o^TFj X! '^gn'^n / Ptor(2fe-i)-i(a;)ptor(2fe)-i(y)sgn(t/-a;)di^i(x)dzyi(y) =Pf Af. 



tESl fe=l' 



4.4.6 Simplifying Bf 

Turning to Bit , we first write 



M 



detVP,,(/3)= Y^ Sgncr Jl Pt,oa(2m-l)-l(/5mHoa(2m)-l(An), 
o"eS2M m=l 



SO that 



(y^S2M m— 1 

1 *^ /• 



<^&S2M m=l' 



Next we set 



n2Af = {o- e S2M ■■ cr(2m - 1) < cr(2TO) for m = 1, 2, . . . , M}, 



then 



1 



M 



^^' " Jjl, Yl Sg^'^n / [Pt'o<T(2m-l)-l(/3K'o<T(2m)(^) -Pt'o<T(2m)(/3K'oa(2m-l)-l(/3)J '^^2(^) 

' o-enaM "1=1 " 

= PfBP. 

4.4.7 The Pfaffian formulation for the total partition function 

It follows from (|4?2T|) that 

ZZ'm"^ E sgnt.PfAf.PfBP, 

i:LyN_ 

and therefore that, 

Z''i'''^(X)= X! ^^ E sgnt-PfAf-PfRP 

= Y. T. sgnt-X^PfAf.PfBP. 

{L,M) i-.L/'N 
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Next, we set X = XI where I is the N x N identity matrix. It is clear that, for each 

(XAPXT)t = XtAfXT, 

and 

Pf(XtAfX^) = detXfPfAf =X^PfAf. 

It follows that 

Z''''"^{X)= Y^ Y^ sgni-Pf(XAPX"^)t -PfBP =Pf(XAPX"^ + BP) 

{L,M) i-.Ly'N 

where the last equation comes from the formula for the Pfaffian of a sum of antisymmetric 
matrices [1^. Finally, since XApX^ — X^A^, we arrive at Theorem 14.31 

4.5 Proof of Theorem ISTSl 

Given xi, . . . , xn, yi, . . . , y^ G M and indcterminants ai, . . . , a^r, bi, . . . ,bN we define the 
measures rji and r]2 on R given by 

N N 

d?7i(a) = ^a„(i(5(a-x„) and d772(/3) = ^ &„d(5(/3 - j/„), 

n— 1 n—1 

where S is the probability measure with unit point mass at 0. Using these measures, we will 
specialize the situation in Section |4^ to the measures i^i = w{^ + rii) and 1^2 = w^ {iJ^ + V^) ■ 
We will derive a Pfaffian form for the correlation functions of the microcanonical ensemble 
with weight w by expanding both the integral and Pfaffian sides of ()4.19p for this choice of 
vi and V2 and equating coefficients of the various products of the indcterminants. 

4.5.1 Expanding the Integral Definition of Z"^'"'^ 
Starting with (|4T7)) . and setting Z"^^""^ == Z''^^''^{1), we have 

(L,M) ' ' "^ " 

It is easily verified that 

L 

d{f^ + Vi)^ia)^Y Y. dr1i{a.^)d^i^-\a^') 
and 

M 

m=0 p:m/'M 

SO that. 



Z 1^ 1^ 1^ 1^ 1^ x!M! 

{L,M) e=Q m=0 u:l^Lo:m^M_ 
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We may relabel the a and the (3 in any manner we see fit. In particular, we may replace 
each u : £ /^ L and : m /^ M_hy i : i /^ L and i : m /^ M_ respectively (the redundancy 
in the notation should cause no confusion). We compensate by a factor of (^)(„) which 
counts the number of pairs (u, o) we are summing over. That is, 

EEE 



(L,M)£=Om=0 V ; V / 



Now, 

IN m N 

dri[{a.i) ^'^'^aud5{aj -Xu) and d-q^ifS^) ^ TT Y^ &t> dd{(3k - Vv), 

J = lU=l fe=lt)=l 

and exchanging the sum and product in each of these expressions, 

£ m 

d?7i(ai) = ^ Y["-u(j)dS{aj ~Xu(j)) and dr^^il^i) ^ ^ Wa„^k)d5{Pk - yx>(k))- 

n-l-^N_j = l v:rn'^N_k=l 

Notice that the sums in the latter two expressions are over all (not only increasing) functions 
from £ and m into N_. 
So far, we have that 

^= T.T.T. T. E £!(L - £)!m!(M - m)! L L' "^■^^"' ^^ 

^ \ n "'uU)d<>iaj - a;u(j)) \dfi^-\cti,)l J| K{k)dSi(3k - y^^k)) Ua*^'^""(A') 

*- j=l ^ ^ k=l ^ 

L M _. f I m - 

= E EE E E lUL - ly.mUM - m)! I n '^"O) IT ^^ 

X / / f7L.M(xuVa,y„V/3)dM''"^(aiOdM'^-'"(A')- 

Next we notice that, if u or are not injective, 

^L,M(xuVa,y„V;3) =0, 

and we can therefore replace the sums over i-^N_ and rn^-N_ with sums over i^-^N_ and 
m^^N_. In fact, since il^^Mi^u V ct, yo V /3) is symmetric in the coordinates of Xu and y„, 
we may replace these sums with sums over £ /^ N_ and m /^ N_so long as we compensate 
by factors of £\ and m\. Putting these observations with the definition of Re.m, we arrive at 
the fact that 

^=J2J2J2 E E n«"0)n^''wh''.™("-^«)- (4.22) 

(L,M) ^=0 "1=0 u:£7'JV Dim/' JV *- j = l fc=l ■' 

That is, Z'^^''^'^ /Z is the generating function for the correlation functions of our microcanon- 
ical ensemble. 



24 



A Solvable Mixed Charge Ensemble on the Line: Global Results 



4.5.2 Expanding the Pfaffian Formulation of Z"^'"^ 
It is easily computed that 



(/l.9)r = (/l.9)i + 2^a„ ^f{xn)ei9{xn) - gixn)eifixn) 

N N , s 
r, V^ V^ J/ \~/ •,Sgn(Xn~Xmj 
-^2^ 2^ anamJ(Xn)g(Xm) ^ , 



and 



AT 

{f\9)T = (/l5)4 + Yl ^" [fiyn)<^2giyn) " giyn)e2fiyn) 



For convenience let us write 



£i,i{x,y) ^ ^sg^iv ~ x) and £i^2{x,y) ^£2,iix,y) =£2,2{x,y) ^0, 



and define 
SO that, 



2a„ 1 < n < TV; 

5„_Ar N <n<2N, 



1 1 < n < TV; 

2 TV < 77 < 2Ar, 

and z„ = 



2N 



Xn l<n<N; 

yn-N N <n<2N. 



(fly)? + {f\g)7 = {f\g)l + (/|5>4 + X! '^" [/(2n)ei(n)5(^") - g{Zn)^i(n)f{Zn) 

n=l 
2N 2N 



n— 1 m— 1 

Defining A^ and B^ as in Theorem 13. II and 

we immediately see that 

AP'"! + BP'''^ = CP +WP, 

=AP+BP 

where the j, k entry of Wp is given by 

2N 
^C„ \pj{Zn)ei(^n)Pk{Zn) " Pk{Zn)<ii(n)P]{Zn)\ 



n=l 



2N 2N 



/ ^ / ^ ^n^'raPj\Zn)Ph\Zrn)^i{ra).i[ri)\Zrn^Z'fi). 



n— 1 m^l 

Next we define the TV x AN matrix X by 
X = 



^mPj \^rn) y ^m^i{7n)Pj \^m) 



j = 0,...iV-l; m = l,...,2A^. 
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We also define the AN x AN matrix J by 






1 






-1 
















1 






-1 






and the iN x AN matrix Y by 



Y= J 







27V 



Finahy, we set Z = (Cp)-t = [QM^k^o- 
A bit of matrix algebra reveals that 



and therefore, 



AP>'^i + B^^"^ = Z-'^ + XYX"^, 



Z"^^"^ Pf(Z-T - XYX"^) 



Z Pf(Z-T) 

This is useful, since by the Pfaffian Cauchy-Binet identity OS], 

Z''u'^2 ^ pf(Y-T _ x^ZX) 



Pf(Y-T) 



A bit more matrix algebra reveals that 



x'zx = 



And since 



Af-l N-1 

j,k=0 j-k=0 

N-1 N-1 

yCnCm /—i ^i{m)Pj\^m)C,j,kPk{Zn) y/CnCm /—i ^i{m)Pj\^m)C,j,k^i(n)Pk(Zn) 
j.k=a j,k=0 



2N 



Y-T = J 







2N 



m,n— 1 



we have that Pf(Y"'^) = (-1)^. This implies that 



Pf fX'^ZX - Y""^) 
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Looking at the entries of X^ZX and using the definitions of >f7v, ci and 62, we see 

AT-l 

y^ Pj{Zm)Cj.kPk{Zn) = XN{Zm,Zn) 
j,k=0 

N-1 

/ , Pj{Zm,)Cj,kii{n)Pk{Zn) = XN^i{n){Zm, Zn) 
].k=0 

/ , <^i(m)PjiZm)C],kPk{Zn) = £i(m)^NiZm, Zn) 



j,k=0 



N-1 



} , <^i(m)Pj{Zm)C],k€-i{n)Pk{Zn) = ei(m)>(N{Zm, Zn)ei{n) ■ 
j,k=0 



It foUows that 



where 



yui ,1^2 



K 



l)'^Pf(-J-K) =Pf(J + K) 



i{m) A{n) , 



-i2N 



/ r r K"- ''"• '(7 7 \ 



Using these definitions and the formula for the Pfaffian of the sum J + K, [TU], we find 



that 



2JV 



E E pfK, 



n=l i:n/^2N 

where Kj is the 2n x 2n antisymmetric matrix given by 



Kt = 



Vct(i)Ct(fe)^)^'^^^''°*^'°'' (^tO) , 2t(fc) ) 



i,fe=i 



For each i: n /^ 2N there exists non- negative integers £ and m such that i? + m = n and 
functions u : £ /^ N_ and D : m /* iV given by 



u(j) = t(j) for aU j such that i(j) < N 



and 



0(.?) = t(i) - iV for ah j such that A^ < t(j) < 2Ar. 

In this situation we write i = uV 0. (We ahow for the possibihty that u or D is the "empty" 
function, in which case uVo = tioruVO = u respectively.) It follows that we may write 



N N 



EE E E pfK„.„ 

e=0 m=0 u:l/^N_0:rn/^N_ 



Finally, 



pf i^uvu = h^ n «"0) n ^"W f pf 

'- j=l k=l ' 



Kji (a;u(j),a;u(jv)) Kj^ ixu{j),yo{k')] 

K^ {yx,{k),Xu(j')) Kj^ {yx>{k),yx>(k')) 



j,j' ,fc,fc' — 1 
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and thus, 



Z 



N N 



EE E E n«uo)n^" 



oik] 



£=0 m=0 u:ey^N_v:my^N_ ^ j=l 

X 2^Pf 



fc=l 



(4.23) 



K^ {yv{k),Xu{j')) K^ (yt.(fc),2/o(fc')). 



j,j',k,k'=l 



Comparing the coefficient of aia2 ■ ■ ■ aibib2 ■ ■ ■ bm in this expression with that in (|4.22p 
we find that 



i?,,,„(x;y)=2^Pf 






j,j',k,k'=l 



as desired. 



(4.24) 



4.6 Proof of Theorem [331 

Let Hn be the standard Hermite polynomial. It is known (cf. [1]) 

{H2m, H2n+1 — 4n_ff2»i-l)l = 4/l2n<5m,n, 

where 

h„ := f Hl{x)e-^^dx = ^Tn\. 

Using the fact that H'^{x) — 2kHk-i{x), it follows readily that 

{H2m,H2n+l)4 = f {H2^{x)H'2,,+ ^{x) - H'2^{x)H2n+l{x)) dx (4.25) 

= 2(2n + l)hn5m,n - 4:mh2m-lSm,n+l 
~ h2n+l5m,n — ^2n+l<J)n,n+l • 

We look for skew orthogonal polynomials in the form of 

rri 

P2m{x) ^'^akH2k{x), ao = 1, 

fe=0 

and determine the coefficients cik, ^ < k < m,hy the orthogonal conditions 

{P2,n,H2k+i - ^kH2k-i) - 0, fc = 0, 1, . . . , m - 1. 
This gives, by (I4.24p and (14. 25^ . the following equations on Uk, 

-h2kak-i + [h2k+i + 4(fc + X'^)h2k] o-k - /^2fe+2afe+l =0, A: = 0, 1, . . . , m - 1, 
where we define a_i = 0. Rescale by setting ak '■= /iVy flfc, the above equations become 
-{k + l)ak+i + (2fc + i - X^)ak - {k + i)Sfe_i, fc = 0, 1, . . . , m - 1, 
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which can be used to determine ak recursively, starting from a_i = and ao = 1. It 

_i 
turns out, however, that this is precisely the three-term recurrence relation for L^ ^ (— X~). 

_i 
Consequently, au = Lj.^ {—X"^). Hence, we conclude 



P2m{x) = ^ akH2k{x), 



ak 



fc=0 



k\ 

W) 



■L-H-X^). 



(4.26) 



The Hermite polynomials are related to the Laguerre polynomials by [111 p. 106] 

H2kix) = {-lf2^''k\Ll^'\x^), H2k+iix) = {-lf2^^+^k\xL~^'\x^). (4.27) 

Using this relation and the facts that 



i^(0) 



k + n\ 9"^^ 

, and (2fc)! = ^r(A:+i)r(fc + l), 



(4.28) 



we see that (|4.26l) becomes p.3p . Since P2m is determined by (|4.24l) and (|4.25p . the same 
process shows that 



P2m+i{x) = ^afe(i?2fc+i(a;) - ^kH2k{x)), 



ak = 



k=0 



k\ 

W) 



■L-H-X^). 



(4.29) 



Using the fact that Hl^Jx) = AjH2j^i{x) and H2j+i{x) = 2xH2j{x) — Hl^Jx), we see that 

rn 
P2rn+l{x) ^^ak{2xH2k{x) " 2i7^Jx)) = 2xP2„^{x) - 2P^„(x). 



fc=0 



It remains to compute {P2m,P2m+i)- By (|4.24p and (|4.25p . we have 

{P2m,P2m+l) = am{P2m, H2m+1 ~ 4m_ff2m-l) 

= am[-h2mam-l + {h2m+l + 4(m + X'^)h2rn)a2m\ 

= a™V^2'"+V! [-(2m - 1)L„^ ^(-X^) + (4m + 2X^ + 1)L™^(-X2) 

= a„v^2™+im!2(m + \)L;X^{-X^), 

where the last step follows from the three-term relation of the Laguerre polynomials, from 
which (13. 5p follows readily. 

Finally, we turn to the proof of (EH). Using £L'^{x) = -L'^+l{x) and L'^{x) = 
Lt\^) - Ll^M (m P- 102]), (E3D gives 

rn 
P2ra+l{x)=2xY,{-l)^LlH-X^) [if (x^) + 4-1 C^^') 

TTl— 1 

Y,{-lfLl{x^)[LlH-X^)-Lll,(-X^))+{-irL^H^X^)Li{x^) 



= 2x 



k=0 



where L^(a;) := i^(a;)/i^(0), the formula dS^]) then follows from 

xLUx)^-{a + l){Lt{x)~Lt^,{x)), 

which is a rescaling of the identity xL'^+'^{x) = -{k + l)L'^_^_^{x) + {k + 1 + a)L'^{x) ([III 
p. 102]). 
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